Impact of tissue microstructure on a model of cardiac electromechanics based on MRI data

Cardiac motion is a vital process as it sustains the pumping of blood in the body. For this reason motion abnormalities are often associated with severe cardiac pathologies. Clinical imaging techniques, such as MRI, are powerful in assessing motion abnormalities but their connection with pathology often remains unknown.

Computational models of cardiac motion, integrating imaging data, would thus be of great help in linking tissue structure (i.e. cells organisation into fibres and sheets) to motion abnormalities and to pathology. Current models, though, are not able yet to correctly predict realistic cardiac motion in the healthy or diseased heart.

Our hypothesis is that a more realistic description of tissue structure within an electromechanical model of the heart, with structural information extracted from data rather than mathematically defined, and a more careful definition of tissue material properties, would better represent the high heterogeneity of cardiac tissue, thus improving the predictive power of the model

Feature detection from echocardiography images using local phase information

Ultrasound images are characterized by their special speckle appearance, low contrast, and low signal-to-noise ratio. It is always challenging to extract important clinical information from these images. An important step before formal analysis is to transform the image to significant features of interest. Intensity based methods do not perform particularly well on ultrasound images. However, it has been previously shown that these images respond well to local phase-based methods which are theoretically intensity-invariant and thus suitable for ultrasound images. We extend the previous local phase-based method to detect features using the local phase computed from monogenic signal which is an isotropic extension of the analytic signal. We apply our method of multiscale feature-asymmetry measurement and local phase-gradient computation to cardiac ultrasound (echocardiography) images for the detection of endocardial, epicardial and myocardial centerline

A poroelastic model coupled to a fluid network with applications in lung modelling

Here we develop a lung ventilation model, based a continuum poroelastic representation of lung parenchyma and a 0D airway tree flow model. For the poroelastic approximation we design and implement a lowest order stabilised finite element method. This component is strongly coupled to the 0D airway tree model. The framework is applied to a realistic lung anatomical model derived from computed tomography data and an artificially generated airway tree to model the conducting airway region. Numerical simulations produce physiologically realistic solutions, and demonstrate the effect of airway constriction and reduced tissue elasticity on ventilation, tissue stress and alveolar pressure distribution. The key advantage of the model is the ability to provide insight into the mutual dependence between ventilation and deformation. This is essential when studying lung diseases, such as chronic obstructive pulmonary disease and pulmonary fibrosis. Thus the model can be used to form a better understanding of integrated lung mechanics in both the healthy and diseased states

Modeling 3D cardiac contraction and relaxation with point cloud deformation networks

Global single-valued biomarkers of cardiac function typically used in clinical practice, such as ejection fraction, provide limited insight on the true 3D cardiac deformation process and hence, limit the understanding of both healthy and pathological cardiac mechanics. In this work, we propose the Point Cloud Deformation Network (PCD-Net) as a novel geometric deep learning approach to model 3D cardiac contraction and relaxation between the extreme ends of the cardiac cycle. It employs the recent advances in point cloud-based deep learning into an encoder-decoder structure, in order to enable efficient multi-scale feature learning directly on multi-class 3D point cloud representations of the cardiac anatomy. We evaluate our approach on a large dataset of over 10,000 cases from the UK Biobank study and find average Chamfer distances between the predicted and ground truth anatomies below the pixel resolution of the underlying image acquisition. Furthermore, we observe similar clinical metrics between predicted and ground truth populations and show that the PCD-Net can successfully capture subpopulation-specific differences between normal subjects and myocardial infarction (MI) patients. We then demonstrate that the learned 3D deformation patterns outperform multiple clinical benchmarks by 13% and 7% in terms of area under the receiver operating characteristic curve for the tasks of prevalent MI detection and incident MI prediction and by 7% in terms of Harrell's concordance index for MI survival analysis

Multi-objective point cloud autoencoders for explainable myocardial infarction prediction

Myocardial infarction (MI) is one of the most common causes of death in the world. Image-based biomarkers commonly used in the clinic, such as ejection fraction, fail to capture more complex patterns in the heart's 3D anatomy and thus limit diagnostic accuracy. In this work, we present the multi-objective point cloud autoencoder as a novel geometric deep learning approach for explainable infarction prediction, based on multi-class 3D point cloud representations of cardiac anatomy and function. Its architecture consists of multiple task-specific branches connected by a low-dimensional latent space to allow for effective multi-objective learning of both reconstruction and MI prediction, while capturing pathology-specific 3D shape information in an interpretable latent space. Furthermore, its hierarchical branch design with point cloud-based deep learning operations enables efficient multi-scale feature learning directly on high-resolution anatomy point clouds. In our experiments on a large UK Biobank dataset, the multi-objective point cloud autoencoder is able to accurately reconstruct multi-temporal 3D shapes with Chamfer distances between predicted and input anatomies below the underlying images' pixel resolution. Our method outperforms multiple machine learning and deep learning benchmarks for the task of incident MI prediction by 19% in terms of Area Under the Receiver Operating Characteristic curve. In addition, its task-specific compact latent space exhibits easily separable control and MI clusters with clinically plausible associations between subject encodings and corresponding 3D shapes, thus demonstrating the explainability of the prediction

Guía para ofrecer feedback a estudiantes

Esta guía está vinculada al proyecto de innovación docente: DISEÑO DE UNA METODOLOGÍA PARA LA COMUNICACIÓN EN EL PROCESO DE EVALUACIÓN DE ESTUDIANTES: EL FEEDBACK PROFESOR/A-ESTUDIANTE (66/FO11/25)El objetivo fundamental de la guía, es ofrecer una metodología muy simple y fácil, para dar feedback dirigido a los estudiantes, donde éstos puedan recabar información de sus profesores acerca de los principales aspectos relativos a la evaluación de sus trabajos, en un plazo razonable de tiempo tras la finalización de un trabajo, actividad o asignatura

Multi-class point cloud completion networks for 3D cardiac anatomy reconstruction from cine magnetic resonance images

Cine magnetic resonance imaging (MRI) is the current gold standard for the assessment of cardiac anatomy and function. However, it typically only acquires a set of two-dimensional (2D) slices of the underlying three-dimensional (3D) anatomy of the heart, thus limiting the understanding and analysis of both healthy and pathological cardiac morphology and physiology. In this paper, we propose a novel fully automatic surface reconstruction pipeline capable of reconstructing multi-class 3D cardiac anatomy meshes from raw cine MRI acquisitions. Its key component is a multi-class point cloud completion network (PCCN) capable of correcting both the sparsity and misalignment issues of the 3D reconstruction task in a unified model. We first evaluate the PCCN on a large synthetic dataset of biventricular anatomies and observe Chamfer distances between reconstructed and gold standard anatomies below or similar to the underlying image resolution for multiple levels of slice misalignment. Furthermore, we find a reduction in reconstruction error compared to a benchmark 3D U-Net by 32% and 24% in terms of Hausdorff distance and mean surface distance, respectively. We then apply the PCCN as part of our automated reconstruction pipeline to 1000 subjects from the UK Biobank study in a cross-domain transfer setting and demonstrate its ability to reconstruct accurate and topologically plausible biventricular heart meshes with clinical metrics comparable to the previous literature. Finally, we investigate the robustness of our proposed approach and observe its capacity to successfully handle multiple common outlier conditions

3D Shape-Based Myocardial Infarction Prediction Using Point Cloud Classification Networks

Myocardial infarction (MI) is one of the most prevalent cardiovascular diseases with associated clinical decision-making typically based on single-valued imaging biomarkers. However, such metrics only approximate the complex 3D structure and physiology of the heart and hence hinder a better understanding and prediction of MI outcomes. In this work, we investigate the utility of complete 3D cardiac shapes in the form of point clouds for an improved detection of MI events. To this end, we propose a fully automatic multi-step pipeline consisting of a 3D cardiac surface reconstruction step followed by a point cloud classification network. Our method utilizes recent advances in geometric deep learning on point clouds to enable direct and efficient multi-scale learning on high-resolution surface models of the cardiac anatomy. We evaluate our approach on 1068 UK Biobank subjects for the tasks of prevalent MI detection and incident MI prediction and find improvements of ~13% and ~5% respectively over clinical benchmarks. Furthermore, we analyze the role of each ventricle and cardiac phase for 3D shape-based MI detection and conduct a visual analysis of the morphological and physiological patterns typically associated with MI outcomes.Comment: Accepted at EMBC 202

MRI-Based Computational Torso/Biventricular Multiscale Models to Investigate the Impact of Anatomical Variability on the ECG QRS Complex

Aims:Patient-to-patient anatomical differences are an important source of variability in the electrocardiogram, and they may compromise the identification of pathological electrophysiological abnormalities. This study aims at quantifying the contribution of variability in ventricular and torso anatomies to differences in QRS complexes of the 12-lead ECG using computer simulations. Methods:A computational pipeline is presented that enables computer simulations using human torso/biventricular anatomically based electrophysiological models from clinically standard magnetic resonance imaging (MRI). The ventricular model includes membrane kinetics represented by the biophysically detailed O’Hara Rudy model modified for tissue heterogeneity and includes fiber orientation based on the Streeter rule. A population of 265 torso/biventricular models was generated by combining ventricular and torso anatomies obtained from clinically standard MRIs, augmented with a statistical shape model of the body. 12-lead ECGs were simulated on the 265 human torso/biventricular electrophysiology models, and QRS morphology,duration and amplitude were quantified in each ECG lead for each of the human torso-biventricular models. Results:QRS morphologies in limb leads are mainly determined by ventricular anatomy,while in the precordial leads, and especially V1 to V4, they are determined by heart position within the torso. Differences in ventricular orientation within the torso can explain morphological variability from monophasic to biphasic QRS complexes. QRS duration ismainly influenced by myocardial volume, while it is hardly affected by the torso anatomyor position. An average increase of 0.12±0.05 ms in QRS duration is obtained for eachcm3of myocardial volume across all the leads while it hardly changed due to changes in torso volume. Conclusion:Computer simulations using populations of human torso/biventricular models based on clinical MRI enable quantification of anatomical causes of variability in the QRS complex of the 12-lead ECG. The human models presented also pave theway toward their use as testbeds in silico clinical trial